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Abstract 

We study direct currents in a simple holographic realization of a superconducting 
film. We investigate how the presence of a DC current affects the superconducting 
phase transition, which becomes first order for any non-vanishing value of the cur- 
rent, as well as several other properties of the superconductor such as the AC conduc- 
tivity. Near the critical temperature we find a quantitative agreement with several 
properties of Ginzburg-Landau superconducting films, for example the squared ratio 
of the maximal and minimal condensate is equal to two thirds. We also comment 
on the extension of our construction to holographic Josephson junctions. 



1 Introduction 



In yy it was shown that a black hole solution in a theory with a charged scalar field 
coupled to Maxwell-Einstein gravity may become classically unstable below some 
critical temperature Tc. This instability induces charged scalar hair for the black 
hole for T < Tc. According to the AdS/CFT duality, the holographic dual of such 
a system is a thermal quantum field theory in fiat Minkowski space with a global 
U{1) symmetry, which is spontaneously broken below by the condensation of 
the operator dual to the bulk complex scalar. In this sense, the boundary theory 
has many of the necessary ingredients to describe a superconductor, or a superfiuid 
[2]. This result was first exploited in [3] to assemble a gravity dual of a system 
undergoing a superconducting phase transition. This construction has been widely 
studied and generalized, and the different gravity duals sharing these same basic 
features go under the name of holographic superconductors (see [U El E] for reviews 
and references). 

In a holographic superconductor, a background magnetic field induces a current. 
However, because of the absence of a dynamical gauge field, this current does not 
expel the magnetic field, unlike the current induced in ordinary superconductors. In 
this sense holographic superconductors more closely resemble thin superconducting 
films or wires. Motivated by this analogy, in the present note we will compare and 
contrast holographic superconductor phenomenology with that of superconducting 
films. More precisely, we study a holographic superconductor in two spatial dimen- 
sions, i.e. a (extremely) thin superconducting film, with a DC current, analyzing 
its phase diagram and the behavior of some interesting thermodynamic quantities. 

In fact, such a system is interesting to study for a second, more ambitious, 
reason. It is believed that holographic superconductors may give an understanding 
of some basic features of high temperature superconductors (HTS). HTS typically 
enjoy a layered structure and, according to the Lawrence-Doniach model p'J, may 
be approximated by films of superconductors separated by Josephson junctions. 
Therefore, a holographic realization of a Josephson junction would be desirable. The 
latter is based on the Josephson effect |8], the phenomenon of current flow across 
two weakly coupled superconductors separated by a very thin insulating barrier (the 
Josephson junction). Our model can then be seen as a flrst, necessary step towards 
the realization of a Josephson junction. 
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We pursue the phenomenological approach of [3] and therefore the gravitational 
system we consider is Einstein-Maxwell theory in four dimensions minimally coupled 
to a charged massive scalar field. In [3] the following basic set up was considered: 
the condensation of the charged scalar in a black hole metric at finite charge density. 
The black hole introduces a temperature T. The finite charge density, which is taken 
care of by allowing for a non trivial profile for the temporal component Aq of the 
gauge field, provides an independent scale needed to get a critical temperature T^. 
One then finds that for temperatures below Tc the charged scalar condenses. Under 
the AdS/CFT map gauge symmetries on the gravity side are dual to global ones on 
the field theory. Then the condensation of the scalar nicely realizes the spontaneous 
breaking of a global U{1) symmetry. 

We want to modify this basic scenario and allow for the presence of a DC current. 
To this end, we consider solutions where also a spatial component of the gauge field 
has a non trivial profile, this providing, via the AdS/CFT map, a current in the dual 
theory. Such solutions can be easily found in the superconducting phase, where the 
scalar is non-zero. However, as pointed out in [9], in the normal state (where the 
symmetry is not broken and the scalar is hence vanishing) the only allowed solutions 
for the spatial components of the gauge field are the trivial ones. For this reason, 
within the minimal Einstein-Maxwell framework, we cannot construct a model de- 
scribing the normal state with DC current. However, as we will discuss in detail, 
this inconvenience will not impede us to obtain some robust results characterizing 
the behavior of holographic superconductors at fixed DC current. 

An additional limitation of our approach comes from the fact that, as in [3], we 
work in the probe approximation. This is the limit where the backreaction of the 
gauge and scalar fields on the metric is neglected. Hence our results are reliable only 
in the regime where the backreaction can be effectively neglected. Luckily, while the 
probe approximation breaks down in the zero temperature limit, for temperatures 
significantly different from zero the results obtained in the probe limit are not sub- 
stantially modified by the backreaction [lOj. There is therefore a large region where 
our results should not be sensibly different from a honest fully back-reacted modeM 

^In [m [12], it has been shown that a phenomenological model of the hkes of [3] can be consis- 
tently embedded in M-theory or type JIB String Theory. Such embeddings constitute an important 
advance towards the understanding of the underlying microscopic theory of the holographic super- 
conductors. Unfortunately, sticking to the probe approximation prevents us from using the models 
of [11] [12] since the charge there is fixed to a finite value, while the probe approximation holds in 
the large charge limit. 
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The same model we are going to study here was already considered in [9l [T3] . 
Differently from those analyses we study the system at fixed current. This choice, 
beside being closer in spirit to real-life experiments, allows us to obtain new results 
about the phase diagram of this system together with interesting checks and predic- 
tions for the behavior of holographic superconductors with DC current, as we now 
summarize: 

• At any finite DC current the transition between the superconducting and the 
normal state is a first order phase transition. We study the temperature depen- 
dence of the condensate and compute the free energy in the superconducting 
state, concluding that at the phase transition the condensate always jumps a 
finite distance to zero. This is a clear indication of a first order phase transi- 
tion. 

• We determine the relation between the current and the superfiuid velocity. 
We largely find nice agreement with expectations for physical superconducting 
films, both for temperatures appreciably lower than, and close to, Tc. More- 
over, it turns out that the form of these curves further justifies the assertion 
in the previous point, namely that the phase transition is first order. Interest- 
ingly, at low temperatures we find that, in contrast with BCS superconducting 
films, for each value of the superfiuid velocity there are two possible values of 
the current. A free energy computation then shows that only one value, in 
fact the highest one, is thermodynamically stable. 

• We study the temperature dependence of the critical current and of the ratio 
given by the value of the condensate at zero current over the value at the 
critical current. For temperatures close to we find that the holographic 
superconductors reproduce the universal results predicted by the Ginzburg- 
Landau (GL) model for superconducting films. On the other hand, at lower 
temperatures our results deviate significantly from the ones of GL. 

• Finally, we study the dependence of the AC conductivity on the DC current. 
We present results for the conductivity in the direction transverse to the cur- 
rent. At low temperatures we analyze the dependence of the frequency gap 
on the DC current. As expected on physical grounds, as we increase the cur- 
rent the frequency gap diminishes, and it does so down to a minimal (but not 
vanishing) value where the first order phase transition occurs. 
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The rest of this note is organized as follows. In section 2 we present the bulk 
Lagrangian, the equations of motion which we have solved numerically, and motivate 
our ansatz and boundary conditions. In section 3 we present the physical output of 
our numerical studies, namely the checks and predictions mentioned above. Section 4 
contains our conclusions as well as a possible strategy for constructing a holographic 
dual of (an array of) Josephson junctions. 



2 The gravity dual of a DC superconductor 

As advertised, we pursue a bottom-up approach to holographic superconductivity, 
and consider as a starting point the model originally presented in [3], Einstein- 
Maxwell theory in 4-dimensions minimally coupled to a charged, massive scalar 
field. We stick to the probe approximation and in this case the action of the scalar- 
Einstein-Maxwell theory reduces to 



S = I dx 



-\f^uF^'' -\{d,- tA^)m\'' - m*m 



(1) 



where the Einstein-Hilbert term has been suppressed, since the backreaction of the 
fields on the metric can be ignored in the probe limit (the Einstein equations de- 
couplejfl. F^i, is the U{1) field strength, ^! is the complex scalar with charge 1 
and mass m, and g is the determinant of the metric g^^, which we take to be the 
asymptotically AdS planar black hole metric 

ds^ = -f\r)dt' + —- + —{dx' + dy^) where /(r) = — - — . (2) 
j[r) r 

The radial direction extends from the black hole horizon at r = tq = (ML^)^/^ to 
the boundary of AdS at r — )■ oo, L is the radius of AdS and M the mass of the 
black hole. Beside the holographic coordinate r, we have three others (t, x, y), which 
parametrize the AdS boundary and hence the (2+l)-dimensional dual field theory 
space-time. 



^ Let us start from the standard Lagrangian \-~jF^i,F'^^ — 1(9^ — i — |\I>P] + 

Einstein, rewrite it in terms of the rescaled fields — q"^ and — qA^ and 

take the limit q — !■ oo while keeping ^I^ , fixed. Then, the Lagrangian becomes 

1 / '' 
^ - 

ter sources decouple from the Einstein equations and the dynamics of the vector and the scalar 



-jF^uFf^" - \{df, - i - j^'p + Einstein. Due to the l/q'^ factor the mat- 
es decouple from the Einstein equations and the dynamics of the vector a 
field are described by the action ([T]) in a vacuum solution of the Einstein equations. 
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The temperature of the black hole (and hence of the dual field theory) is given 

by 

As in [5] we will take the scalar mass to be = —2/ LP' which is above the 
Breitenlohner-Friedman bound. 



2.1 The ansatz 

According to the AdS/CFT map, the VEV of the U{1) current in the dual field 
theory is identified with the subleading boundary asymptotics of the bulk gauge 
field. Hence, to describe holographically a superconductor with a DC current, we 
need to look for bulk solutions where the black hole develops charged scalar hair 
in the presence of a non-trivial profile for a spatial component of the gauge field. 
More precisely, we are interested in getting a current in the x direction, therefore we 
will look for solutions which are independent of the time coordinate t and of y, but 
with a non-trivial dependence on both r and x. We choose the gauge = (this 
leaves the freedom to perform r-independent gauge transformations, as we will do 
later). As our solutions are y- independent, we set Ay = 0. Thus we must determine 
A^, At, and \E' as functions of x and r. 

We choose the modulus of the scalar to be independent of x and similarly for A^ 
and At. However, having a current then requires that the phase of be x-dependent. 
The simplest such ansatz reads 

^(r,a;) =^(r)e^^% (4) 

which automatically satisfies the equations of motion for Ar and for the phase of 
^E', the latter imposing that 6 is indeed a constant. Summarizing, tp, A^ and At 
are functions of r while 6' is a constant. Notice that in a superconductor the spatial 
derivative of the phase of the condensate is the superfiuid velocity [HI [13]. Therefore, 
in our case we are describing a superconductor with a constant superfiuid velocity, 
since by the AdS/CFT map the latter is given by 6. 

Notice that a non-zero Ar^ contributes positively to the effective scalar mass, 
hence one expects that a sufficiently large A^ will win against the negative contribu- 
tion coming from the time-component of the gauge potential, eventually destroying 
black hole superconductivity This corresponds to a critical maximal current 
in the dual field theory, above which the system enters the normal phase, which is 
indeed what is expected for physical superconductors. 
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The equations of motion following from the action ([T]) are 

dr{r^drAt)-2'-j^At = 0, (5) 

dr{fdrA,)-2^' {A, -6) = , (6) 

dr{r'fdr^)-L^A,-ef^+'^^ + ^^ = 0. (7) 



These equations are invariant under two independent scaling symmetries 

r Xr , ro -> Aro , L L , ip ^ ip , At ^ XAt , A^ XA^ , 9 ^ X9 , (8) 
r ^ r , ro -)■ ro , L ^ uL , u^'^ip , At -)■ v'^At , A^ -> i'~'^A^ , ^ z/"^^ . 

When performing numeric computations we find it convenient to work with variables 
and coordinates which are invariant with respect to these rescalings (and hence 
dimensionless) . In our case they are 

- , —A, , — , LV'. (9) 

ro ro ro 

The equations of motion written in terms of these rescaled and dimensionless quan- 
tities are the same as before with ro and L set equal to one. 

2.2 Asymptotics and their dual interpretation 

The equations of motion ([5])- ([7]) are second order, and so we expect six constants 
of integration, which together with 9 imply seven parameters. Regularity at the 
horizon sets = at r = 1 and this, via the equations of motion, imposes two 
more constraints on ip and A^. This leaves four parameters. In addition, there are 
boundary conditions at the boundary of AdS. The leading asymptotics of the fields 
at large r read 

^. = Af-^+0(r-^) , At = Af^-^+0{r-^) , ^ = ^ + ^+0(r-3) , 

(10) 

while 9 is constant everywhere. The leading contribution of the time and space 
components of the gauge field correspond, via the AdS/CFT map, to a chemical 
potential n and a source for the x-component of the dual current, respectively. This 
source for the current will be nothing else than the superfiuid velocity Vx- Since the 
gauge field is covariantly coupled to the scalar, on the boundary we get a term of the 
form dx ^ — A^x"* =9 — A^x'^ and we see that one can, through a gauge transformation. 
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trade 9 for Ax IT3] . In fact, from now on we will choose to work in the gauge 
^ = and identify the superfluid velocity with Af\ The subleading asymptotics 
correspond instead to the charge density p and the VEV of the current density J^. 

With our choice of scalar mass term we have two options for the corresponding 
dual operator. This is because both asymptotic behaviors of the scalar are nor- 
malizable at the boundary, so both of them can correspond to a VEV of a dual 
operator [15]. We can choose the leading asymptotic coefficient ^^^^ to be a source 
of a scaling dimension 2 operator O2. In this case the expectation value of O2 will 
be proportional to the subleading coefficient ■ip^'^\ Or, we can choose ■ip^'^^ to be a 
source of an operator Oi with scaling dimension 1. In this case the VEV of Oi 
will be proportional to For definiteness, in what follows we will work with the 
operator O2 (the basic results are not qualitatively different for the opposite choice). 
To have spontaneous breaking of the U{1) symmetry, we want the source to vanish 
and hence we will impose ip^^^ = 0. 

In summary, we have seven parameters, three regularity conditions at the horizon 
plus the two conditions ip^^^ = , ^ = 0. Hence, we expect a two-parameter family 
of solutions, which may be parametrized, for example, by the temperature and 
current of the superconducting film. Given the temperature and current one may 
then calculate the value of the order parameter (O2). 

Undoing the rescaling described above, we can rewrite the field theory quantities 
in terms of the asymptotic coefficients of the dimensionless fields. We find 

,= fTA?\ p={J,)='-^T^A^\ 

= f T , ^ ( J.) = A« , (O2) = 2^ , 

(11) 

where we have used the field and coordinate redefinitions given by eq. (Q, and 
written tq in terms of the temperature via eq. ([3]). 

As can be seen from the expressions for the chemical potential and charge density 
in eq. f llip . the asymptotic behavior of At only determines the dimensionless ratios 
fi/T and p/T^- In other words, the gravity dual only gives us information about the 
dimensionless ratio of the two scales in the theory: the temperature and the charge 
density (or the chemical potential). One can then decide to use either the chemical 
potential or the charge density to fix a scale. We will use the former and study the 
evolution of {O2) / p^ as a function of T/p and jxlp^- Accordingly, using eq. ( fTTj) . 
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we define Tc, the critical temperature at zero current, as 

^ = -^ (12) 

where A^°''|c is the critical value of Af'^ for which the condensate turns on at jx = 
0. When studying the thermodynamics of the system we will work in the grand 
canonical ensemble, which corresponds to a system at fixed chemical potential. 



3 Holographic predictions 

We have numerically solved the system of coupled differential equations ©-(IT]) and 
determined the condensate (O2) as a function of the current and the temperature. 
Using a shooting technique we have integrated the equations from the horizon up 
to the boundary, with the boundary conditions discussed before. Via the AdS/CFT 
maps detailed in eqs. ffTTj) . we have then determined the surface of solutions for the 
condensate as a function of the current and the temperature. 



3.1 Phase transition with DC current 

The first important thing we want to analyze is how the presence of the current 
modifies the temperature dependence of the condensate. This is shown in figure 
[1] for different values of the current and compared with the result at zero current 
obtained in [3]. Two significant modifications occur. First, at any finite value 
of the current one can see that the curve (O2) vs T becomes bivaluated. There 
appears a new branch (the dotted line) corresponding to states where the value of 
the condensate is much lower. In the following, by computing the free energy we will 
see that the states with lower value of the condensate have a larger free energy than 
their counterparts with larger (O2) at the same temperature. Therefore, this new 
branch corresponds to thermodynamically disfavored states B Second and more 
importantly, we observe that the superconducting state exists up to a maximum 
value of the temperature (where the plot turns back). Crucially, at that point the 
value of the condensate is larger than zero. Therefore, at the phase transition the 
condensate must jump a finite distance to zero. Unless one fine tunes the parameters, 
such a jump will almost certainly change the energy and so require some latent heat, 

^Notice that this branch smoothly joins the (unstable) normal phase branch of ^ in the Jx ^ 
limit. 
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Figure 1: On the left we plot {O2) versus the temperature for several values of the 
current: from the innermost to the outermost jx/T'^ = 28.98, 14.49, 2.90, 0.290. The 
dotted lines correspond to the states with larger free energy than their counterparts at 
the same temperature. On the right we show for comparison the result at zero current. 
Notice that at the critical temperature (O2) vanishes in this case. 

implying that the phase transition is first order. Moreover, as is expected on physical 
grounds, the temperature at which the phase transition occurs is always lower than 
Tc, the critical temperature at zero current, and its value decreases with increasing 
current. 

This phase transition pattern is quite different from that of refs. [21 [13] • The 
analysis performed there corresponds to experiments where instead of the current, 
the superfluid velocity is kept fixed. There it was found that the superconducting 
phase is separated from the normal phase by a second order phase transition from 
zero superfluid velocity up to a tricritical point where the phase transition becomes 
first order and remains so up to the maximum velocity, where the phase transition 
would be at zero temperature (similar results were found in [16], in the context of 
superconducting D-brane models). In fact, as we will see in section [321 the different 
phase transition pattern one finds when working at finite current agrees with what 
is known about the relation between the current jx and the superfluid velocity in 
superconducting films. 

3.1.1 The free energy 

In order to confirm our previous claim, namely that the states with lower value 
of the condensate are metastable, we shall now compute the free energy of the 
superconducting phase and show that it is larger for the metastable branch (dotted 
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line in figure [T]) . 

The free energy of tlie system is determined by tlie action ([1]) evaluated on-shell 
Q = —T Sos plus possible boundary counterterms [IT] . For the present case the reg- 
ularized action was presented in [13] . We shall proceed along those lines to compute 
the free energy of the physical configuration we are interested in. Substituting the 
equations of motion ©-(IT]) into the action one finds 



^0 



(fx 



+ I d'x[i^'Al-j^'A 



(13) 

which is the unregularized on-shell action (the prime means derivative with respect 
to r). This action consists of three boundary terms resulting from the kinetic terms of 
the temporal and spatial components of the gauge field, and the scalar, respectively; 
plus a bulk contribution coming from the interaction terms. From the asymptotic 
behavior of the fields ( ITOl) it follows that only the boundary term corresponding to 
the scalar field if) is divergent, and thus we need to add the corresponding coun- 
terterm. Moreover, one must specify the boundary conditions which are imposed at 
infinity on the various fields. In our case, one should add boundary terms which take 
us to an ensemble where v4|°^ and A^x^ are held fixed, corresponding via eq. ( ITTl) 
to (O*-^-*), the chemical potential jj, and the current jx- AH in alL the boundary term 
that does the whole job reads (see [H] for a rigorous analysis )□ 



(fx {r^ + 2r^ tfjtfj' + Ax A'^) 



(14) 



Substituting the behavior of the fields written in eq. ffTOj) into the regularized on- 
shell action given by the sum of the contributions ( IT3l) and ( |T4l) yields the following 
expression for the free energy 



\{aT Af^ + Af A^^) - ^^^^^ ^^-^ 



where V stands for the volume of the system. 

^Varying in the bulk yields a boundary term ^ A'^SA^, which provides a good variational 
principle for a boundary condition (5A^|r=oo = 0. This corresponds to an ensemble where we 
are keeping fixed the asymptotic value of A^, which in our case means fixed chemical potential 
At\r=oo) and fixed source of the current (~ Ax\r=oo)- One can go to an ensemble where 
5S ^ A^5 A'^ by adding a boundary term ~ r'^A^ AJ^|r=oo- 
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Figure 2: For a value of the current jxl^c — 2-9, we plot the free energy of the supercon- 
ducting phase zooming in on the region in which T is closer to its maximum value (where 
the plot turns back). The dashed line corresponds to the points with a lower value of the 
condensate at a given temperature. We show on the right the corresponding plot of (O^^^) 
versus T. One can see that the lower branch (dotted line) corresponds indeed to states 
with larger free energy and thus metastable. 

We can now compute the free energy of the superconducting states making up 
the plot in figure [1] and confirm that the lower branch (dotted line) is metastable. 
This is shown in figure [2l Notice, however, that we cannot determine precisely at 
which value of the temperature the phase transition occurs. In order to know this, 
we would need to compare the free energy of the superconducting state with the free 
energy of the normal state at the same value of the current. Unfortunately, within 
the minimal framework we are using, it does not seem possible to describe such a 
normal state. Let us elaborate a bit more on this. Naively, the first thing one could 
try to do is to look for a solution with non-trivial At and A^^. but vanishing scalar 

= 0). However, as noticed in [9], the only such solution satisfying regularity 
conditions at the horizon has = identically. This result should be expected on 
physical grounds. In the normal state the dual system is no longer superconducting 
and thus one expects that in the absence of an electric field the current must vanish 
(the DC conductivity is now finite). One could then try to switch on a background 
electric field in the x direction through the addition of a contribution of the form 
—E ■ t to A,j.. However, in the Maxwell action having Ft^ = ~E does not modify 
the equation of motion for A^ which is then still divergent. One would expect 
that using a richer model describing non-linear interactions between the bulk fields 
could solve the problem. Indeed, as shown in [19] and more recently in [20], a DBI 
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action does provide a solution with non-vanishing current in the normal state. There 
the conductivity was computed and found to depend both on the electric field and 
the charge density. Yet it is not clear how to implement the scalar condensation 
corresponding to the superconducting phase in this scenario A different approach, 
closer in spirit to the model we are dealing with here, consists in going beyond the 
probe approximation, thus looking for a solution which in the normal state would 
correspond to an asymptotically AdS charged black hole with vector hair (with At{r) 
and = -E ■ t + h{r)). 

Let us emphasize that although within the probe approximation regime we can- 
not compute the free energy of a normal state with current, the conclusion about the 
phase transition being first order is robust. It is clear from our computations (figure 
[1]) that at the maximum temperature the value of the condensate is different from 
zero and then it must jump during the phase transition. If we were able to compute 
the free energy of the normal state it may be that the phase transition would occur 
at a value of the temperature somewhat lower than the maximum value. However, 
as we see from the plot, the condensate would still be different from zero at that 
point. Simply, the superconducting state would be metastable from the actual tem- 
perature of the phase transition up to the maximum temperature (as it happens for 
instance in [9| 113]). 



3.2 Current and velocity 

In this section we study the relation between the current jx and the superfluid 
velocity v^- As we explained at the beginning of section [3], the integration of the 
equations ([5])- ([7]) results in a two-parameter family of solutions, which we chose to 
parametrize in terms of T and jx- This means that once T and jx are fixed all other 
physical quantities of interest are determined up to a discrete choice, in particular 
also the superfiuid velocity Vx. We will now fix T and obtain a one-parameter curve 
of solutions relating jx and z/^.. The result is presented in figure [3] for several values 
of T. Close to the critical temperature (left panel) the relation is an upside down 
paraboloid which becomes smaller as the temperature approaches Tc (eventually 
shrinking to a point for T = Tc). On the other hand, for low temperatures (right 
panel) the relation between jx and Vx is linear almost all the way up to a given 
maximum velocity above which the superconducting state exists no more. 

In general our results match nicely with what is known about the relation be- 

^For a p-wave superconductor such a DBI construction is possible, see for instance |21) . 
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Figure 3: Plots of the current jx versus the superfluid velocity I'x at fixed temperature. On 
the left panel we show the results for two temperatures close to the critical temperature: 
T = 0.998 Tc (solid line) and T = 0.996 Tc (dashed line). On the right we present the 
curve we find for T = 0.244 Tc, the dashed line corresponds to metastable states since they 
are solutions with higher free energy than their counterparts with the same value of the 
current or the velocity. 

tween the current and the superfluid velocity in thin superconducting films [H]. As 
we comment below, though, one qualitative difference with respect to BCS super- 
conducting films is that at low temperatures at the maximum superfluid velocity 
the current is non-zero. 

For thin films at temperatures close to T^, where the GL model is reliable, the 
jx versus Ux curve has exactly the same features as the one in the left panel of 
figure El As we will now explain this is responsible for the different phase transition 
pattern one finds when working at fixed current or at fixed superfluid velocity, 
respectively. For every value of the current there are two values of the superfluid 
velocity. The current is clearly zero at zero velocity, but also at the maximum value 
of the superfluid velocity the current falls to zero because the condensate vanishes at 
this velocity (this holds also in the present holographic model, see [9l|13]). Therefore, 
at the maximum value of the current the superfluid velocity is not at its maximum 
and the condensate has a finite non-zero value. This means that if one increases 
the current a bit further there is no corresponding value of the superfluid velocity 
in the superconducting phase. Then, the system passes into the normal phase with 
the condensate jumping a finite distance to zero, which, as already explained, is a 
signal of a first order phase transition. At the critical temperature this argument 
breaks down, as no current is possible in the superconducting phase and there is no 
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discontinuity. There the phase transition is second order. 

At low temperatures the relation between and for thin films is linear from 
zero up to a maximum velocity, the depearing velocity, at which the current falls 
steeply to zero [H]. As we see on the right plot of figure [3] we find this linear behavior 
for a long range of currents. However, unlike BCS films, we also find that at the 
maximum value of the superfiuid velocity the current is non-vanishing. According 
to the analysis in the previous section, for that value of the current the condensate 
is non-zero and hence if one increases the superfiuid velocity the condensate will 
jump to zero. Hence the phase transition is first order. This agrees with the result 
of P [13] where it is found that the phase transition at low temperatures is indeed 
first order. 

As already noticed, for each value of the current one finds two solutions with 
different values of the velocity and so two values of the condensate (O2). As 
shown in figure [2], we have found that the free energy calculated from the gravity 
solution using eq. (fT5|) is always lower for the configuration in which the magnitude 
of (O2) is higher. 

3.3 Critical current and critical condensate 

In this section we will discuss two further results of our holographic analysis. Recall 
that our model aims to describe a thin superconducting film and that it is expected 
to be reliable for the whole range of temperatures (except for very low temperatures, 
where the backreaction needs to be taken into account). For temperatures near Tc 
the GL theory is expected to give an accurate description of such physical system 
and therefore we have to compare with GL in this regime. It is important in the 
GL derivation that the film is thin, as this allows one to ignore the free energy 
contribution of the magnetic field generated by the current. As our magnetic field is 
non-dynamical, it will not be generated by a current, and so it will not contribute to 
the free energy. This makes the comparison between holographic superconductors, 
which are inherently ungauged, with GL model particularly sound for thin films. 
On the other hand, for temperatures far below the critical temperature our results 
give new insights on the phase diagram of holographic superconductors. 

GL theory predicts that near Tc the critical current jc is proportional to (Tc — 
T)^/^. As illustrated in figure H] we find that this scaling is indeed obeyed by holo- 
graphic superconductors for temperatures close to Tc. On the other hand, at low 
temperatures our results differ appreciably from GL scaling. This is to be expected. 
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Figure 4: Plot of the critical current versus the temperature. The left panel shows a 
log-log plot from which we can read-off the critical exponent, getting 1.497, which agrees 
with the expected GL scaling of 3/2 within our numerical precision. The right panel shows 
the departure from GL scaling (solid line) at low temperatures. 

For large currents the temperature at which the phase transition occurs is apprecia- 
bly lower than T^, and hence far from the regime where the GL effective description 
is valid. Moreover, being the phase transition first order, one does not expect any 
power-law scaling. Notice that we are defining the critical current to be the highest 
current at which the superconducting solution exists, at a given temperature. As 
we discussed in section I3.1.H it might be that the phase transition happens for a 
lower value of the current and thus the value we are considering would correspond 
to a metastable state. Nevertheless, at temperatures close to Tc one can reasonably 
expect than the corrections to the free energy coming from the very small value 
of the current are almost negligible and the critical current agrees with the max- 
imum allowed value. The fact that under this assumption the result got for the 
holographic superconductors agrees with GL is an a posteriori reassurance. Being 
more conservative, one should take our result as an upper bound, especially for the 
large current (corresponding to low temperature) regime: in other words, at a given 
temperature, the corresponding critical current would be at most equal to the one 
predicted by the plot in figure HI 

A second prediction of the GL theory is that, at any fixed temperature, the 
norm of the condensate monotonically decreases with respect to the velocity from 
its maximum value (02)00- The critical current is reached before the maximum 
velocity, when the norm of the condensate has an intermediate value {02)0- More 
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Figure 5: Plot of the ratio {{02)c/ (02)00)^ versus tlie temperature. Tlie solid line 
corresponds to the value of 2/3 predicted by the GL theory. 



precisely one has 

We have found numerically that this relation is also satisfied for holographic super- 
conductors at temperatures near the critical temperature Tg. This can be seen in 
figure O where we plot the ratio f lTB]) versus the temperature. Again, away from 
the behavior changes sensibly. Notice that the same warning about our inability to 
determine exactly the critical current applies here and therefore, especially in the 
region T ^ Tc, the points in figure |5] should be considered as a lower bound for the 
ratio {{02)0 1 (02)00)^- Notice that the ratio goes to a constant at zero temperature. 



3.4 Conductivity 

In this section we will study the AC conductivity of the system and characterize its 
dependence on the DC current. To compute the conductivity one must consider an 
electromagnetic perturbation on top of the hairy black hole solution. This is easy 
for a perturbation along the direction orthogonal to the current [i.e. a perturbation 
of Ay), since it decouples from other perturbations of the gauge vector or the scalar 
field. Conversely, a perturbation of couples to perturbations of A^ and i\} and 
hence the computation of the conductivity along the direction parallel to the current 
becomes more involved and we will not attempt to do that here. 
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The equation of motion for a zero-momentum perturbation 6 Ay = e *'^*y4y(r) 
takes the form 



where we have apphed again the rescahngs given in eq. and defined cD = 
3/{A7c) u/T. Notice that this is the same equation as considered in [3J, but the 
background solution for the scalar ip is different now, in particular it depends on the 
current. The boundary asymptotics (r — > oo) of Ay take the form 



The conductivity is given by the zero-momentum retarded current- current correlator 
which by using the AdS/CFT dictionary can be calculated in terms of solutions 
satisfying ingoing wave boundary conditions at the horizon [22]. In fact we recover 
Ohm's law on the boundary 



where we have taken into account that Ay is introducing a background potential 
on the boundary and thus an electric field Ey = —dt Af\ 

By solving numerically eq. ( |T7I) with infalling boundary conditions at the horizon 
(r = 1) we can compute the conductivity as a function of the frequency u at given 
values of temperature and current. In figure [6] we show the results obtained at a low 
temperature (T = 0.04 ■ Tc) for different values of the current. The AC conductivity 
displays the features already observed in [5]. At large frequencies it approaches a 
constant, a characteristic of theories with AdS^ duals [26]. On the other hand, at 
a; = we expect a delta function in Re(o"), a fact confirmed, via the Kramers- 
Kronig relations, through the presence of a pole in the imaginary part of o" at a; = 0. 
Finally, we see that for small enough frequencies, within our numerical precision, 
Re((T) vanishes. This gap can be parametrized in terms of a critical frequency Ug. 
As we show in figure [6] there is a minimum of lm{a) around the point where Re((T) 
becomes non-zero. Then, following |23], we define Ug as the frequency minimizing 
the imaginary part of the conductivity. 

For weakly coupled superconductors the gap is predicted to be oUg/Tc = 3.5 at 
T = f]3]. In the original model of [3] the gap was found to be oJg/Tc ~ 8 and it 
seems that such a high value holds quite generically for holographic superconductors. 
This has been seen as an indication that holographic superconductors are indeed 




(17) 






(19) 
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Figure 6: On the left we plot the real part of the conductivity versus the frequency for sev- 
eral values of the current at T = 0.04 T^. From left to right: j^/T^ = 43.6 , 41.3 , 29.6 , 2 • 
10~^. The leftmost curve corresponds to the maximum current at this temperature. On 
the right we show both the real (solid line) and imaginary (dashed) part of the conductiv- 
ity as a function of the frequency again at T = 0.04 Tc and with jx/T"^ = 43.6, the largest 
allowed current at that temperature. 



strongly coupled. In figure [7] we plot OJg as a function of the current. In the region 
of very low current we recover the result Wg/Tc ~ 8. As we increase the current, 
ujg decreases continuously until we reach the maximum current where it has a finite 
value. This is consistent with the phase transition being first order at that point. 
The condensate is non-vanishing and thus we expect Ug to be also different from 
zero. 

In weakly coupled superconductors a definite relation exists between Ug and the 
energy gap A at zero temperature, Ug = 2 A, A(T) being the minimum energy 
required for charged excitations at a given temperature T. In strongly coupled 
superconductors one does not expect the gap to necessarily satisfy this relation so one 
could have wondered what the relation is for holographic superconductors. However, 
as noticed in [101 EH ES] and recently reviewed in [6] , holographic superconductors 
are not hard-gapped, in general, and a non-zero conductivity is present even at 
small frequencies, though exponentially suppressed (for recent work on hard-gapped 
holographic superconductors see for instance [27], and [28] for a recent discussion 
on this point). Indeed, computations of the temperature dependence of the specific 
heat [21] showed that holographic superconductors behave similarly to some strongly 
coupled superconductors as heavy fermion compounds: the specific heat does not 
vanish exponentially at low temperature (this being a consequence of, and hence an 
indication for, the existence of an energy gap), but as a power law. 
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Figure 7: Plot of LOg as a function of the current. At zero current we find w^/Tc = 8.87, 
while at the maximum current i^g/Tc = 6.28. 

4 Conclusions and outlook 

In this paper we have considered a simple holographic model of a thin superconduct- 
ing film with DC current. We focused on the modifications that the presence of a DC 
current induces on the thermodynamics as compared to the same holographic model 
with vanishing current, originally studied in [3]. Most notably, the phase transition 
becomes first order for any finite value of the current. Moreover, the conductivity 
gap becomes a function of the current, too: the frequency gap diminishes as one 
increases the current but never reaches zero before the phase transition occurs, in 
agreement with the phase transition being first order. Other results we obtain nicely 
agree with expectations for thin superconducting films, such as the relation between 
the current and the superfiuid velocity, both at low and high (that is near to Tc) 
temperatures. The only qualitative difference is that at sufficiently low tempera- 
tures and high superfiuid velocities, the velocity no longer uniquely determines the 
current. 

Ideally one would like to go beyond the probe approximation. Besides leading 
to better control over the very low temperature regime, this is necessary in order 
to obtain a holographic description of the normal phase. To describe the supercon- 
ductor in the normal phase with a DC current one should switch on an external 
electric field, which is needed to keep a constant current in the normal phase, since 
there we have a non-vanishing resistivity. However, as we already noticed, only the 
full system of coupled scalar Maxwell-Einstein gravity equations could in principle 
allow for a non-trivial but meaningful (that is non-singular) solution. This implies 
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that one should go for a fully back-reacted analysis. 

4.1 Towards holographic Josephson junctions 

Our primary motivation for this work was the observation that HTS's have typically 
a layered structure and may, according to the highly successful Lawrence- Doniach 
model [7], be approximated by superconducting films separated by Josephson junc- 
tions. While the present model can be seen as a first step in this direction, one 
would like to find a complete holographic description of a Josephson junction. Let 
us elaborate a bit on this. 

One kind of Josephson junction that appears particularly amenable to a holo- 
graphic construction is the S-c-S junction. Such a junction is composed entirely 
of the same superconducting material, but the superconductor is thinner at the 
junction, and so for example will have a lower critical current. Imposing a space- 
dependent metric in the boundary theory, one might be able to cook-up a dual model 
with varying critical current. This might seem quite ad hoc in a bottom-up context 
but it might possibly arise quite naturally in a string theory context. For example, 
one may consider M-theory compactified on a 7-dimensional Sasaki-Einstein man- 
ifold fibered over AdS^. The critical current, like all quantities in a holographic 
superconductor, depends on the compactification geometry. Then, to construct a 
junction, one may merely need to vary the moduli of the 7-manifold in the region 
corresponding to the junction, that is, over a finite interval in one of the field the- 
ory directions. An array of junctions would then correspond to moduli that vary 
periodically in one field theory direction. 

Ideally such a spatial dependence of the moduli will be a solution of the super- 
gravity equations of motion. One hope of realizing such a solution is as follows. 
Each Josephson function may correspond to a brane extending from the boundary 
to the horizon and along the field theory directions parallel to the Josephson junc- 
tion, and also potentially wrapping some cycle of the compactification manifold. 
One may then hope that by correctly choosing this cycle, one may engineer a ge- 
ometry in which the backreaction of the brane causes the desired deformation of 
the compactification manifold, reducing the critical current and therefore forming 
an S-c-S junction. A Lawrence-Doniach HTS may then correspond to an array of 
such branes. 
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